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Abstract: In this paper we compute the holographic entanglement entropy for massive fla¬ 
vors in the D3-D7 system, for arbitrary mass and various entangling region geometries. We 
show that the universal terms in the entanglement entropy exactly match those computed 
in the dual theory using conformal perturbation theory. We derive holographically the 
universal terms in the entanglement entropy for a CFT perturbed by a relevant operator, 
up to second order in the coupling; our results are valid for any entangling region geometry. 
We present a new method for computing the entanglement entropy of any top-down brane 
probe system using Kaluza-Klein holography and illustrate our results with massive flavors 
at finite density. Finally we discuss the differential entropy for brane probe systems, empha¬ 
sising that the differential entropy captures only the effective lower-dimensional Einstein 
metric rather than the ten-dimensional geometry. 
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1. Introduction 

In recent years there has been considerable interest in entanglement entropy and its holo¬ 
graphic realisation, following the proposal of [1] that entanglement entropy can be com¬ 
puted from the area of a bulk minimal surface cohomologous to a boundary entangling 
region. This proposal was proved for spherical entangling regions in [2] and arguments 
supporting the Ryu-Takayanagi prescription based on generalised entropy were given in 
[3]. Entanglement entropy has been computed in a wide range of holographic systems, see 
the review [4]. 

The focus of this paper is on the computation of holographic entanglement entropy in 
top-down brane probe systems, which are widely used in phenomenological applications of 
holography. Entanglement entropy is a new computable for such systems and, following 
the pioneering works of [5, 6], can act as an order parameter for confinement and other 
phase transitions. 

A top-down brane probe system is expressed in terms of a ten-dimensional supergravity 
background and a brane embedding into this background. The Ryu-Takayanagi prescrip¬ 
tion is however based on extremal surfaces in the reduced Einstein {d + l)-dimensional 
metric, where d is the dimension of the dual field theory. One of the main results of 
this paper is a systematic method to compute the holographic entanglement entropy for 
any top-down brane probe system, using the method of Kaluza-Klein holography [7] to 
extract the lower-dimensional Einstein metric. This method reproduces earlier results of 
[8, 9, 10, 11] but allows entanglement entropy to be computed for any brane system with 
arbitrary worldvolume fluxes. (Earlier results for massless flavors at finite density can be 
found in [12].) We illustrate our methodology using the example of the D3-D7 system at 
finite mass and density. We compute the holographic entanglement entropy for massive 
flavors, with arbitrary mass and various entangling region geometries, and use our new 
methodology to address the case of finite density. 

Brane systems provide a new testing ground for the dependence of entanglement entropy 
on the field theory and on the shape of entangling region, topics of considerable current 
interest, see for example [13, 14, 15, 16, 17, 18]. In particular, one can explore the structure 
of universal logarithmic terms; these are well-understood for conformal field theories (see 
e.g. [19, 20]) and recent papers have explored the behaviour of entanglement entropy under 
relevant perturbations using conformal perturbation theory [21, 22, 23, 24, 25, 26]. It was 
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shown in [23] that for a CFT deformed by a relevant operator 




d'^xO 


( 1 . 1 ) 


there is a new logarithmic divergence in the entanglement entropy of the half space 


5S 


M\^ 


(d — 2) TT 

4(d- 1)1^ 


d+2 

2 


d+2 

2 


^log 



( 1 . 2 ) 


when A = (d + 2)/2 with A the area of the dividing surface and N the normalisation of the 
two point function of the operator O. Here euv and em correspond to UV and IR cutoffs 
respectively. 

In section 6 we prove (1.2) by analysing the volume divergences of the holographic en¬ 
tanglement entropy and show that (as postulated in [23]) such a divergence occurs for an 
entangling surface with arbitrary geometry. We also show explicitly that (1.2) agrees with 
the logarithmic terms in the entanglement entropy for the D3-D7 system at finite mass, 
using the results of [27] to determine the holographic two point function normalisation. As 
well as matching the universal terms in the entanglement entropy, we explain the origin of 
finite terms in the entanglement entropy for massive flavor systems, in terms of the effective 
IR description of the system in terms of a CFT deformed by irrelevant operators. 


There is a growing literature connecting quantum entanglement with the global structure of 
the bulk spacetime, see in particular [28, 29]. In [30] a relation between the area of generic 
(non-minimal) surfaces and entanglement was proposed and this idea was sharpened with 
the introduction of differential entropy [31, 32, 33, 34, 35, 36]. We verify that the differential 
entropy in the D3-D7 system indeed computes the area of a hole in the reduced Einstein 
metric; the agreement is somewhat subtle since the depth of the hole is itself corrected by 
the presence of the probe branes. 

In section 8 we discuss the implications of the fact that the entanglement and differential 
entropy are related to the reduced Einstein metric, rather than the ten-dimensional metric: 
even if entanglement allows us to reconstruct the reduced Einstein metric completely, this 
information does not suffice to reconstruct the ten-dimensional geometry. Moreover, the 
causal structure in ten dimensions only agrees with that of the reduced Einstein metric in 
special cases (e.g. product metrics); the global structure is qualitatively different between 
five and ten dimensions even for well-understood examples such as the Coulomb branch of 
N = A SYM. Reconstruction of the full ten-dimensional geometry would therefore seem to 
require a generalized notion of entanglement in the dual field theory. 


The plan of the paper is as follows. In section 2 we briefly review the relevant features 
of the D3-D7 system. In sections 3 and 4 we compute the entanglement entropy for slab, 
half space and spherical entangling regions for the massive D3-D7 system. In section 5 we 
present a general method to compute the entanglement entropy in any brane probe system 
using Kaluza-Klein holography and illustrate our method with the D3-D7 system at finite 
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mass and density. We discuss the field theory interpretation of our results in section 6 and 
give a holographic proof of (1.2) for generic entangling regions. In section 7 we show that 
the differential entropy computes the area of a hole in the Einstein metric and we discuss 
the meaning of entanglement and differential entropy for top-down solutions in section 8, 
illustrating our discussions with Coulomb branch geometries. We conclude in section 9 and 
various technical results relevant to section 5 are contained in appendix A. 

2. Massive flavors 

In this paper we will explore entanglement entropy for massive brane systems, focussing 
for the most part on the specific example of the D3-D7 brane system. Consider W D3- 
branes and Nf <C W parallel coincident D7-branes. As discussed in the early days of the 
AdS/CFT correspondence [37, 38] the decoupling limit gives rise to AA = 4 SYM coupled to 
Nf massless flavors; the resulting field theory is an AA = 2 SCFT. Taking the background 
AdS ^ X metric to be of unit radius: 

ds'^ = ^ -b dx^dx^) -b dO'^ + sin^ OdQ^ -b cos^ Odcj)^ (2.1) 

the embedding of a probe D7-brane corresponding to a massless flavor is described by cj) 
constant and 9 = 7r/2, i.e. the probe D7-brane wraps AdS^ x S^. 

Suppose one separates the D7-branes from the stack of D3-branes; the resulting open 
strings are massive and the field theory in the decoupling limit corresponds to AA = 4 SYM 
coupled to Nf massive flavors [39] (we discuss the massive deformation of the M = 2 SCFT 
further in section 6). The corresponding D7-brane embedding in AdS^ x is described 
by (p being constant and the angle 6 depending on the radial coordinate z as 

sin^d = (I-m^z^), (2.2) 

where m corresponds to the flavor mass, or equivalently the separation of the D7 and D3 
branes. Note that the probe brane wraps the equator of the 5^ as z —)• 0 and smoothly 
caps off at a finite value of z = 1/m, controlled by the flavor mass. 

The D3-D7 brane system can be used to model mesons holographically. Considerable work 
has been done on generalizations of the probe brane embeddings to finite temperature and 
finite density, see the review [40] and on the backreaction of the flavor branes onto the 
geometry [41, 42, 43]. In particular, note that interesting meson melting phase transitions 
are observed at finite temperature and density, see for example [44, 45]. Backreacting mas¬ 
sive flavors is non-trivial even at zero temperature and density, since the flavors break the 
global symmetry to 50(4) and the resulting ten-dimensional metric is therefore of coho¬ 
mogeneity three. Smearing the branes over the compact space reduces the cohomogeneity 
of the metric but this is obscure from the field theory perspective as it corresponds to an 
averaging over different field theories. 

In this work we will calculate the entanglement entropy and the differential entropy for 
the massive flavor system at zero temperature and zero density, and match our results 
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Figure 1: Illustration of a slab boundary region: Ax = I is the width of the region, and is the 
regularized area of its boundary faces. 


with field theory results based on conformal perturbation theory. We will also present 
a method to compute the entanglement entropy for any probe brane system (with or 
without worldvolume gauge fields) and illustrate this method with the case of massive 
flavors at finite density. The method is equally applicable at finite temperature, although 
at finite temperature the entanglement entropy will include both thermal and quantum 
contribution; matching with field theory results is considerably harder as few results for 
finite temperature exist. It would however be interesting to explore the finite temperature 
results in the context of melting phase transitions. 


3. Entanglement entropy for slabs 


In this section we compute the entanglement entropy for a slab on the boundary (see Figure 
1), working to leading order in the ratio of the number of flavors to colors, Nf/Nc- 

Let us begin by reviewing the computation of entanglement entropy for a slab in AdS^. 
We define a slab region on the boundary of width Ax = Z by x G [0,1], and take as an 
embedding ansatz z = z{x) on a f = 0 hypersurface. Defining the regularised lengths of the 
other spatial directions as L, it is then easy to show that the Ryu-Takayanagi entanglement 
entropy functional [1] for the embedding surface is; 


S 


L2 ^ y/IT^ 
4Gn Jo 


(3.1) 


where Gyr is the Newton constant. Since we chose the AdS^ to have unit radius, the 
Newton constant is dimensionless and can be related to the number of colors W as 


1 _ 
SttGn 


(3.2) 
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Figure 2: The minimal surface for a slab boundary region - the boundary is at z = 0 and z = z is 
the turning point of the surface. 


The Lagrangian is independent of x explicitly and hence the associated Hamiltonian is a 
constant of motion. Rearranging the expression for this constant of motion one easily finds: 


z' = 


\/z^ — 


(3.3) 


where 5 is clearly the turning point of the solution since z'{z') = 0. The entanglement 
entropy is then obtained by substituting this solution into the entropy functional, resulting 
in 


2;3 


S = 


L^z 


dz 


N 




(3.4) 


Here we have included a factor of two, from the two halves of the entangling surface, i.e. 
0 < X < //2 and 1/2 < x < 1. It is useful to define the dimensionless parameter s = z/z so 
that z G [e, z] —)• s G [a, 1] where a = e/z is also dimensionless. We obtain, for example: 


dx 

dz 


Vl - s® 


Note that the entanglement entropy is thereby manifestly dimensionless 


S = 




ds 





(3.5) 


(3.6) 





V^r(-1) \ 

6r(i)z2 ) ’ 

where in the latter equation we retain only terms which are finite or divergent as the cutoff 
e ^ 0 . 


It is simple to find the induced metric of the entangling surface and its associated stress 
tensor. The induced metric is given by 7 ^™ = daX^^di,X’''g^,^ where a,b = (s, y, w) run 
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over the surface indices and //, v run over all AdSs indices. The induced metric is therefore 


7aT = 


1 1 


s2(l — 56^ ’ g2-^2 1 g2^2 


We note for further use that 






^/l - 


Differentiating the action functional, the stress tensor for the surface is given by: 


_ _ 2 

— - 


which evaluates to; 


TZn = 


(3.7) 

(3.8) 

(3.9) 

(3.10) 


after using the chain rule. It is then a simple matter to calculate these components, 
resulting in: 

TZn={s\l-s^)z\^,s^z\sH\sH^) (3.11) 

which we will make use of below. Note also that the relation between the width of the 
slab, I, and the turning point of the minimal surface is 


I = 2 


z^dz 


Jo (56-26)2 

3.1 Flavor contribution 


3 - = 25 


fi s^ds 20fr(|) 
/o(i-.6)|- r(i) 


z. 


(3.12) 


Now let us compute the change in the entanglement entropy caused by the presence of iVj 
flavors branes, with Nf <C W- A priori, to compute this change one would expect that 
one needs to compute the backreaction of the branes to linear order in Nf/Nc and then 
extract from the backreacted geometry the change in the five-dimensional Einstein metric 
and hence the change in the area of the minimal surface. As mentioned earlier, it is hard 
to compute the backreacted ten-dimensional geometry because of the high cohomogeneity 
of the problem; smeared solutions are known and entanglement entropy was computed for 
these smeared solutions in [46]. 

It is important to note however that entanglement entropy is defined in terms of the hve- 
dimensional Einstein metric, not the ten-dimensional (Einstein-frame) metric. One does 
not in general obtain the correct answer for the entanglement entropy by computing the 
area of a minimal surface in the ten-dimensional metric, see section 5. Note that [46] used 
the ten-dimensional metric rather than the five-dimensional Einstein metric. 

Computing the full ten-dimensional backreaction without smearing and extracting the 
effective five-dimensional Einstein metric is intractable for general probe brane systems. 
Several methods have therefore been developed to extract the entanglement entropy from 
the probe brane embedding, see [ 8 , 9, 10, 11]. The methods of [ 8 , 10, 11] are particularly 
applicable to spherical entangling regions, for which the CHM map [2] may be exploited. 
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In this section our discussion will follow that of [9], which is applicable to all entangling 
region geometries. 


On general grounds the change in the entanglement entropy for any perturbation in the 
five-dimensional Einstein metric is 

where is the perturbation in the five-dimensional Einstein metric, is the energy 
momentum tensor of the minimal (entangling) surface in the background and the integral is 
over the original entangling surface. Therefore one can compute the entanglement entropy 
provided one can extract the change in the five-dimensional Einstein metric. For general 
brane embeddings the computation of the perturbation in the five-dimensional Einstein 
metric is subtle; in section 5 we present a method to compute the Einstein metric for all 
types of brane embeddings. 


It was observed in [9] that the perturbation in the five-dimensional Einstein metric is 
straightforward to compute whenever the brane embedding has an induced worldvolume 
metric which is diagonal (a product of a non-compact part and a compact part which is 
embedded in the sphere part of the background geometry) and the non-compact part of 
the metric has no dependence on the sphere coordinates. In such a case the linearised 
backreaction on the metric for probe branes can be computed [10] as: 

-Hz),h{z), h{z), h{z)) (3.14) 


where the metric perturbation is sourced by the effective brane energy momentum tensor 


^efF 


i.e. 




(3.15) 


This effective stress energy tensor is obtained by reducing the brane action over the three- 
sphere: 


I =—Tf f (f’a f = —Tj ( (i^cr(27r^)(l — m^z^) 2 (3.16) 

JAdSz Js^ JAdSs 

where 7 q,^ is the worldvolume metric for the brane^. The worldvolume metric is diagonal 
for the given embedding and therefore the determinant factorises, allowing the integral over 
the three-sphere to be evaluated. The resulting effective action then depends only on the 
non-compact part of the worldvolume metric but note that the effective tension of 

this brane is z dependent. Varying this effective action with respect to the non-compact 
part of the background metric results in the effective energy momentum tensor 




(3.17) 


Note that this method for computing the effective source term for the five-dimensional 
Einstein metric relies on the fact that the worldvolume brane metric is a direct product of 

^Note that we denote the worldvolume metric for the brane as 7 and the induced metric on the entangling 
surface as 



Ax 

i -► 



Figure 3: The relationship between the original minimal surface and the D7-probe embedding - 
when z < 1/m the entire minimal surface lies within the probe brane embedding. 


non-compact and compact parts. The method is also not applicable for brane embeddings 
in which worldvolume gauge fields are non-zero or worldvolume fields source other super¬ 
gravity fields as well as the metric. In section 5 we will discuss a more generally applicable 
method for computing the entanglement entropy contributions from probe branes which 
does not rely on a diagonal worldvolume metric. 

Substituting (3.17) into (3.15) gives the following equation: 

fiz) + zh'iz) = fiz) = ^(1 - (3.18) 

Here to is the backreaction parameter, proportional to the number of flavors A//: 

to = WttGnTo; To = (3.19) 

where Tj is the tension of a D7-brane. Only the gauge invariant combination f{z) is 
determined by the Einstein equations. However, continuity of the metric and of the extrinsic 
curvature at z = 1/m requires that h{z) satisfies h{l/m) = /i'(l/m) = 0 . 

Substituting the metric perturbation and the minimal surface stress energy tensor into 
(3.13) thus gives 

J dsdwdy^^^-^^^=== {f{zs){l - s^) -b h{zs){s^ + 2 )) (3.20) 

where this integral is over the original entangling surface with coordinates {s,w,y). Defin¬ 
ing a = L'^/{AGnz"^) for convenience (where we have computed the trivial y, tc-integrals to 
give the factor and taken into account the factor of two arising from the two halves of 
the entangling surface), and unpacking f{zs) we find: 

6S = a [ ds( ~ f{zs) - zh'izs) + (3.21) 

J \ s^^/l - s6 J 
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Integrating the second term by parts, the bulk contribution cancels the third term and one 
is left with a boundary contribution: 


5S = a 




Rzs) 


h{zs) 




(3.22) 


where a = ejz and 6 = l,l//rfor/i< 1,^> 1 respectively. Here ^ = mz is the dimension- 
sless mass parameter. This latter distinction occurs because, although the integral runs 
over the original entangling surface which has 2 : G the integral in fact only receives 

a non-zero contribution when 7 ^ 0, i.e. for 2 ; < 1/m. When 2 < 1/m (i.e. ^ < 1), 
the entire entangling surface lies within the brane embedding, and the upper limit is thus 
z = z or s = 1. When 2 > 1/m (i.e. ^ > 1) however, the upper limit will depend on the 
mass and be given by 2 = 1/m i.e. s = 1/fj, (see Figure 3). 

For both cases, the boundary term at s = b actually vanishes. The expression within square 
brackets trivially vanishes at s = 1 for the case < 1 , and it vanishes for /x > 1 using the 
continuity condition h{l/m) = 0. The expression for both cases is thus given by: 

6S = ai ds ^ ^3 f{zs) + h{e y I (3.23) 


with b depending on the case as mentioned above. Expanding for / this becomes 


5S 



(1 - + h{e) 


VT^ 


(3.24) 


Note that the entanglement entropy depends explicitly on the gauge fixing for the metric 
perturbation. One choice of scheme would be to set f{z) = 0, corresponding to Fefferman- 
Graham coordinates while a second natural choice of scheme is to fix h{z) such that the 
cutoff is unchanged to linear order and one then obtains the relation: 


For this gauge choice one obtains 


6S 


tpL'^ 

A8G]\f 




(i-/iV)2 + 




A8G]\f 




(l-/xV)^ 




(3.26) 


where we note that the gauge dependent contribution bSgauge is independent of the turning 
point 2 , since h(0) is finite. We will discuss this point further below. In what follows we 
will retain the gauge dependence explicitly, rather than fixing a gauge, and show that this 
gauge dependence drops out of universal terms. 
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In computing the integral we first specialise to the case of small mass so b = 1. Performing 
the integral directly over the range s G [a, 1] and expanding the answer in a gives the 
following up to 0(a}: 


ds 


_ ,.2 2\2 _ I ^,,2 I I ,,4 r(l/3) 

^ > 2a2^3^ 12 r(7/6) ^ 


12 r(ll/6) 3 


- log2 + 2^2loga 


(3.27) 


The result thus becomes: 


toL"^ 


SS = 


2 , r(-i/3) 


1 2 
_ _ 1 _ _ TD^ -I- 

48GAr V2e2 3 12^2 r(7/6) 

“l“^5*gauge(^? ^)* 




(3.28) 


We next consider the case of large mass so b = l/fJ.- The result is given in terms of 
generalised hypergeometric functions: 

/"I/a* ../T7TT6 I 1 

J = —+ —3F2({1/2,1,1},{2,2},1/^‘^) (3.29) 

(-1/2,-1/3, 2/3,1/^6) 

+ ^2 Ti (-1/2,1/3,4/3,1/^6) + 2fi^\og{iJa) + 0{a^) 
Expanding for large mass one then obtains 

“ = + 2m=log(m£) - + O (p) + O (^) ^ £). 

(3.30) 

Note that the power and log-divergent terms agree for /U < 1 and ^ > 1. 


We can immediately obtain the change in the entanglement entropy for the half space 
from the 5 —)■ oo limit of the above expression. In this case the entangling surface extends 
throughout the bulk and has no turning point. The contribution to the entanglement 
entropy from the brane is then 

Note that the divergent terms differ from (3.30) by an overall factor of two, since the 
entangling surface in the field theory no longer has two disconnected parts. We will discuss 
the field theory computation of (3.31) in section 6. 


3.2 Changes in turning point and entanglement surface 

The perturbed entangling surface has a turning point for which the relation between the 
turning point and the width of the slab Ax = I is changed relative to (3.12). The equation 
for the perturbed entangling surface is obtained analogously to (3.3) and given by 

=(5 ■ 0 (^ “ 0 “ (5 ■') ■ 


11 - 



I 



' ~ -X- ' ' 
Zt 


Figure 4: Illustration of the change in the turning point of the minimal surface - 5 is the turning 
point of the original minimal surface (i.e. the one that is actually used to compute the flavor 
contribution to the entanglement entropy (3.13)), whereas Zt is the turning point of the minimal 
surface in the backreacted geometry. 


for some constant z, and where f{z) and h{z) are the metric perturbations discussed 
previously. The width of the slab is then given by 

r z^dz ^ w- _ 1 

2 Jo (z® — 2:^)2 2 Jq (^ 6 _^ 6^2 2 

where the upper limits of integration are explained as follows: for /r > 1, the surface ends 
in the region in which the perturbations f{z) and h{z) vanish; from (3.32) the turning 
point therefore remains at z = z, and we continue to use the definition /r = mz. 


z^h 


(z® - z®)2 


(3.33) 


The boundary term in (3.33) vanishes since h{l/m) = 0 and h{0) is finite and therefore the 
relation between the turning point and the slab width depends only on the gauge invariant 
combination of metric perturbations. Substituting this combination using (3.18) one hence 
obtains 


where note that [i 
in 


2 ^ 


I ni) 


+ 




(1 - S ®)2 


dss' 


(3.34) 


mz depends on z implicitly. The integral can be computed resulting 



dss^ 


(1 — S®)2 


1 

2 V 


^16^^(-/r^ + V^® - 1) + 6 2 F 1 


/I 2 5 1 \ 
V2’3’3Vy 


(3.35) 


and expanding for /r ;§> 1 one finds 


l_^J VnTil) to \ 

2 ^ y r(i) 576m^z^ j ' 


(3.36) 


For the small mass case the situation is more complicated, since from equation (3.32) we 
find that the turning point of the surface is itself changed (see Figure 4). Let the perturbed 
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turning point be 


zt = z + to5z. 


(3.37) 


The latter is computed by setting z' = 0 in (3.32), resulting in 

t^dz = -zh{z). (3.38) 

Note that the shift in the turning point depends on the metric perturbation h{z) explicitly, 
rather than the gauge independent combination. 

The relation between I and zt is now calculated using the relation (3.32) which can be 
rewritten as 


= ((1 + hCz))izf - Z^) - fiz)iz! - z^) + H{z){iz^ - z^)y , (3.39) 


where we define 

h{z) = H (z) + h{z), (3.40) 

with H{z) by construction being zero. In relation (3.39) we implicitly work to first order 
in tO) which in particular implies that z can be replaced by zt in terms multiplying f{z) 
and H{z), which are already of order to- Therefore 

1 1 r ^^dz ^ 1 r\3^ Af+^H') 1 

2 {l + h{z))^Jo {zf — z^)^ 2 Jo (zf —2;®)5 2 

The boundary term vanishes at z = 0 and the contribution at z = zt is zero since H{zt) = 0 
to order to- Since h'{z) = H'{z) the combination appearing in the second integral is the 
gauge invariant combination as before and therefore 


z^H 


(^6 _ ^ 6)2 


(3.41) 


l_ 

2 


zt{l - \Hz)) 


s^ds Zt 

(1 _s6)| ^ 24^°, 


ds- 


(l-s6)l 


where now fx = mzt- Computing the integrals we obtain 


/ 


zt{l-\h{~z))^-^^^ 


+ 


tpzt f \/^r(|) _ V 


12 


V r(i) 


+ 


6 r(^) )' 


(3.42) 


(3.43) 


Note that h{z) = h{zt) at this order. 

Even though we only needed the original entangling surface to compute the entanglement 
entropy above, these changes to the turning point are important to keep track of when 
comparing the differential entropy to the gravitational entropy of the corresponding hole 
in the bulk, as we will discuss in section 7. 
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3.3 Finite contributions 


To understand the infra-red behaviour of the entanglement entropy it is often useful to 
isolate the finite contributions. 


A dimensionless, cut-off independent quantity was defined in [47, 48] by differentiating with 
respect to the mass as 

= (3.44) 

note that this expression is valid for four-dimensional quantum held theories, with different 
expressions being proposed in lower dimensions. Implicitly, e and z (or equivalently 1) are 
held hxed. A priori, since the gauge dependent terms depend on the mass it is not obvious 
that this quantity will be independent of the gauge. However, on general grounds the gauge 
dependent terms must make the form 


1^2 

'^S’gauge = CL—2—2 —h CLqW ?-|- (!I(e^) 


(3.45) 


with a _2 and oq dimensionless coefficients. This form follows from the fact that the en¬ 
tanglement entropy is extensive, i.e. it is proportional to the area of each slab L?', and 
the underlying theory is conformal. This implies that the only hnite terms in the scheme 
dependent part of the entanglement entropy must be proportional to m^L^, since m is the 
only other cutoff independent scale in the problem. Since neither a _2 nor oq contribute to 
(3.44), the quantity computed by (3.44) is indeed independent of the gauge. 


Computing this quantity one finds that for < 1 


5Sm. = 


( \/^r(i 


while for ^>1 


SSm. = 


48G]\f y 6r(^ 

toL"^ 


3 / 4 ~2 

—mz 


32Gn 3^2 

which can be expanded for 1 

6S - ^ 


2Fi[-1/2, 1/3,4/3,1/^1 


48Giv y 


— 


1 


8m*z^ 


+ 0 




(3.46) 


(3.47) 


(3.48) 


For slab geometries an alternative method of defining a cut-off independent quantity is (see 
for example [49, 50, 51, 52]^) 

oc 

= (3.49) 

This quantity is manifestly independent of the coordinate choice h{z), since (JS'gauge is 
a local quantity, which is hence independent of the (non-local) slab width I, as we saw 
explicitly below (3.26). For /j, < 1 this quantity evaluates to 


_ toL^ ( ^ 

A8Gn \ 6z2 r(|) ^ 6r(^) "" ^ 

^Si is always positive and decreasing in two dimensions and plays the role of a c-function. 


(3.50) 
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while for 1 one obtains 


48Gn (miosis)) ' 

The limit ^ ^ 1 probes the IR of the theory: for fixed m this corresponds to taking an 
entangling surface which extends deep into the bulk. Therefore a finite quantity should in 
this limit decouple from UV physics. Comparing (3.48) and (3.51), the first quantity does 
not fulfil this criterion (as the term of order m? derives from the logarithmic divergence) 
whereas the latter quantity does. We will hence use (3.51) in section 6 when discussing the 
IR physics. 


3.4 Phase transitions 


The flavor contributions to the entanglement entropy (3.28) and (3.29) match at mz = 1, 
i.e. when the turning point of the entangling surface is at the location where the metric 
perturbation vanishes. This matching was guaranteed by the continuity of the metric 
perturbation and its derivative at z = 1/m. There are however discontinuities in the 
derivatives of the flavor contribution at mz = 1, induced by the discontinuities of higher 
derivatives of the metric perturbation at z = 1/m. In particular, there is a discontinuity in 
the fourth derivative of the entanglement entropy with respect to the slab width (at fixed 
mass): 



(3.52) 


From the field theory perspective it is more natural to fix the mass (i.e. the theory) and 
vary the slab width (i.e. the entangling region). However, if one instead looks at the 
variation of the entanglement entropy with respect to the mass at fixed slab width, the 
fourth derivative is also discontinuous: 


d^S \ 


(3.53) 


Correspondingly the finite quantities SSm and SSi have discontinuities in their second and 
third derivatives, respectively. 

The discontinuity arises from the discontinuity in second derivatives of the metric at z = 
1/m. The fourth derivative of the entanglement entropy contains the terms: 


1 


liV — 






1/m 


(3.54) 


Note that the discontinuity does not arise at lower order in derivatives since the volume 
element of the entangling surface vanishes at the turning point. 


At first sight one might try to assign a physical interpretation to the discontinuity of the 
entanglement entropy, i.e. a phase transition. However, the discontinuity is inherited from 
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the discontinuity in the metric derivatives and correspondingly in the curvature. This 
discontinuity is likely to be an artefact of the probe approximation: in a fully back-reacted 
solution for the D3 and D7 branes there should be no source terms in the energy momentum 
tensor and hence no discontinuities in the curvature of the metric. In other words, one 
would expect from gauge/gravity duality that the backreacted solution should solve the 
type IIB equations with no sources. The metric and curvature should hence be continuous 
and the metric for the back-reacted solution should be smoothened around z = 1/m, over 
a radial coordinate range Az <C 1/m. 

It is interesting to note that at any finite density the brane probe extends throughout the 
bulk and therefore there is no longer any discontinuity at finite z: as we show in section 
5 the backreacted solution is indeed smoothened around z = 1/m, over a small but finite 
radial coordinate range. 


4. Entanglement entropy for spherical regions 


In this section we compute the entanglement entropy for the case of a spherical entangling 
surface, extending the small mass results of [10] to generic mass. The methods of [8, 10, 11] 
are in principle applicable to spherical entangling regions but in practice the CHM map 
[2] becomes intractable for finite mass, as the probe brane embedding in the hyperbolic 
black hole is extremely complicated. Therefore again our discussion will follow closely the 
method of [9]. 


Writing the boundary metric in spherical coordinates we have: 

ds^ = + dr^ + r‘^d^li +dz^) (4.1) 

for the AdS^ metric. We define a ball on the boundary by r < ii and take as an embedding 
ansatz z = z{r) at t = 0. The functional for the entangling surface is then: 


S = 


TT 


v/r 

dr - 


y/2 


GnJo z^ ■ ^ 

where we have done the trivial integral over the two-sphere. It is easy to show that the 
resulting equations of motion are solved by the hemisphere and the desired 

extremal surface is thus a hemisphere in r and z of radius R that wraps the 2-sphere S'^ - 
this surface will be parametrised by {s, 112} where s is defined hy z = Rs and r = R\/l — 
and 112 = {6, </)• 


We can now compute the induced metric on this extremal surface using = daX^^di)X'^ 
where now a,b = (s, 6, (p). One finds: 

^min _ ( _}_ ~ \ /A o\ 

lab -(^g 2 (l_ 52 )’ ^2 (4.3) 

where we note for further use that: 

/l ^ 

(4.4) 
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We can then compute the stress tensor of the surface and one finds (computing only the 
diagonal components since this quantity will be contracted with which is diagonal): 


TZn={s^R\^-S 




1 — 1 — 


rcosec 




(4.5) 


We are now in a position to compute the entanglement entropy, but we must first write 
the metric backreaction in the coordinate system (z,t,r,9,(p): 

= ^^'^^^U{z),h{z),h{z),h{z)gs2) (4.6) 

The resulting entanglement entropy becomes 

((«' + ‘^MRs) - {s^ - l)f{Rs)) (4.7) 


where a = ejR and b = l,l//i for fj, < l,ij, > 1 respectively as for the slab, where now 
fi = mR. Expanding out f{Rs) and again using partial integration on the h\Rs) term as 
in the slab case one obtains: 


-55 = / ds^ — ^P^f{Rs) - "" 


2Gm 


2Gn 


h{Rs 


(1-s2)3/2 


-lb 


(4.8) 


again reducing the contribution of h{z) to a boundary term. The term at s = 6 vanishes 
for both possible values of b for the same reasons as before, leading to: 


-55 = 


tpiT 

2AGn 




3/2 


(4.9) 


Since h{e) depends upon the gauge choice, we can rewrite this expression as in the previous 
section as 

—(1 - (ms)^)^+ -55gauge(e,^,m). (4.10) 

The gauge dependent contribution depends in this case on all three parameters: the cutoff e, 
the mass m and the radius of the spherical region R (note that in the previous expression the 
mass dependence is contained implicitly in the metric function h{z)). Note the difference 
relative to the case of the slab: since the dual theory is local, the gauge dependent terms 
for the slab cannot depend on the slab width. The radius of the sphere however relates to 
the intrinsic curvature of the entangling region, which is a local quantity and therefore can 
appear in the gauge dependent terms. In particular, since h{0) is finite, the non-vanishing 
terms in -55gauge will be either quadratic in R or independent of R in any scheme. 

Let us now consider the small and large mass cases separately. For /U < 1 the contribution 
from the s = b limit to the integral vanishes where 6 = 1. One therefore obtains 


-55 


toTT (R^ 4:^? 4-3 e 


1 + ^ + 

4 9 15^ 


T -55gauge(C) Rj m). 


(4.11) 
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Using the same regularisation scheme as before to fix h{e) one obtains as in [10]: 


6S 


8Gjv(^2e2^ 6 °^2i? 4^ 3 ^ 15 j ’ 


(4.12) 


where 


'^'S'gauge (C) R} 


TTto / R? 2 2 3 

24Gjv ~ 3^ “ 2 ^ 


(4.13) 


For ^ > 1 the extra contribution from the s = b limit of integration is given by; 


tpTr 

24Gjv 



+ o(3 + 4/r )log ( n + 1-^ 




(4.14) 


This vanishes at ;U = 1 as we would expect by continuity. Thus the total contribution to 
the entanglement entropy for ^ > 1 is given by: 


<55 = 


toVT I R^ 4//2 + 3 e 1 8^2 m4 


(4.15) 




8 83/i2 


3 V /u2 \^15 

+(55gauge (^^) Rj TTi). 


+ 


30 


+ 




+ g(3 + 4:^^)l0g ^ ~ ^ 


For fj,^ 1 this expression asymptotes to 


(55 = 


tp-JT 

8Gn 


1 


i?2 /2^2 1\ 3 

+ 2) + 8 - 48^'^ 




T <55gauge(C) Rj 'rn,')' 

(4.16) 

As for the slab, the expressions for the entanglement entropy match at /x = 1, i.e. when the 
turning point of the entangling surface reaches z = 1/m. Derivatives of the entanglement 
entropy with respect to R at fixed m or with respect to m at fixed R become discontin¬ 
uous at /U = 1 because of the metric discontinuity. For spherical entangling surfaces the 
discontinuity arises at fifth order i.e. 


d^S 


m=llR 


(4.17) 


is discontinuous. The discontinuity again arises from the discontinuity in second derivatives 
of the metric at z = 1/m, and is hence expected to be absent in a fully back-reacted solution 
without sources. The fifth derivative of the entanglement entropy contains the terms: 

^ (4.18) 

Note that the discontinuity does not arise at lower order in derivatives since the terms 
contracted with the metric perturbation are zero when the turning point lies at z = 1/m, 
and their first derivative is also zero; see the form of the integrand in (4.8). 
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4.1 Finite contributions 


For a spherical region one can define a finite quantity by differentiating with respect to the 
mass, (3.44). For // < 1 this gives 


For 1 one obtains 

5Sm 


TTto f ^ 
8Gn V 3 4 



(4.20) 


The quantity (3.49) is not finite for a spherical region and it is proposed to use instead 
[53, 54] 


Slm = 

oR 



S. 


(4.21) 


Note that this quantity vanishes for all terms which are independent of R or quadratic in 
R, which in particular guarantees that the gauge dependent terms drop out of Slm- For 
/i < 1 one obtains 


SSlm 


8G]\f 



(4.22) 


For 1 one obtains 


SSlm = .j!"" 2 + 
48GAr^^ 


with all terms with higher order powers in ^ cancelling. 


(4.23) 


As for the slab the limit 1 probes the IR of the theory: for fixed m this corresponds to 
taking an entangling surface which extends deep into the bulk. Therefore a finite quantity 
should in this limit decouple from UV physics. Comparing (4.20) and (4.23), the first 
quantity again does not fulfil this criterion (as the term of order m? derives from the 
logarithmic divergence) whereas the latter quantity does. We will hence use (4.23) in 
section 6 when discussing the IR physics. 


5. Entanglement entropy from Kaluza-Klein holography 

In this section we describe a new method for computing the entanglement entropy of probe 
brane systems using Kaluza-Klein holography [7]. This method is applicable to any probe 
brane system, i.e. for any shape entangling region with any worldvolume gauge fields, and 
can also be used for other systems such as Coulomb branch geometries. 

The holographic entanglement entropy for any static asymptotically anti-de Sitter geom¬ 
etry is given by the Ryu-Takayanagi functional in terms of the area of a minimal surface 
in the Einstein frame metric. Probe brane systems are however described by top down 
constructions. In other words, we first specify a ten-dimensional supergravity solution for 
which a holographic interpretation is known, the usual examples being geometries which 
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are asymptotic to the products of anti-de Sitter and Sasaki-Einstein manifolds. The probe 
system is then specified by the brane embedding into the ten-dimensional background and 
the worldvolume fields on the brane. The backreaction onto the ten-dimensional supergrav¬ 
ity solution is computed by viewing the D-brane action as sourcing the supergravity fields, 
with the sources being localised on the brane embedding. Computation of the backreaction 
therefore involves solving all of the ten-dimensional supergravity equations. 

Even after computing the backreacted ten-dimensional supergravity solution, one can¬ 
not immediately compute the entanglement entropy, because the latter requires the five¬ 
dimensional Einstein metric. Eor any supergravity solution which can be viewed as a 
perturbation of anti-de Sitter cross a Sasaki-Einstein manifold the method of Kaluza-Klein 
holography can however be used to extract the five-dimensional Einstein metric [7]. This 
method implies that, if one only wishes to compute the entanglement entropy, it is not 
actually necessary to compute all of the backreaction of the brane onto the ten-dimersional 
supergravity helds; one only needs to know the backreaction for those field components 
which contribute to the five-dimensional Einstein metric. 

In the rest of the section we will describe the computation of the entanglement entropy using 
the Kaluza-Klein holography approach for massive D7-branes in an AdS^ x background. 
At the end of the section we will discuss further applications and generalisations of this 
method. 

5.1 Kaluza-Klein holography 

The backreaction of the D7-branes onto AdS^ x results in a supergravity background 
which can be expressed as a perturbation of AdS^ x S^. Thus the metric can be expressed 
as 

ds‘^ = idMN + hMN)dx^ dx^ ; (5.1) 

= ^{dz^ + dx^dx^) + {dO'^ + sin^ 6d^l + cos^ 9d(j)^) + hMN{x^, 6a)dx^dx^, 

where we denote ten-dimensional indices as x^] 0 “ collectively denote the five sphere 
coordinates and x™ denote the hve-dimensional coordinates, i.e. {z,x^). Thus g^N 
background AdS^ x metric and Hmn is the metric perturbation. The other type IIB 
supergravity fields are the dilaton (j), the NS-NS three form field strength Hmnp and the 
RR field strengths Fm, Fmnp and Fmnpqr- Only the self-dual five-form field strength 
has a background profile: 

Fmnpqr = Fmnpqr + Imnpqr; (5.2) 

1 , 

= —^dz A dt A dw Adx Ady + sim 9 cos 9d9 A dQ^ A d(j), 
z^ 

with /mnpqr being the perturbation of the hve form field strength. Our normalisation 
conventions are that the Einstein equations for type IIB supergravity are given by 

Rmn = qFmpqrsFj^^^^ H- (5.3) 
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The Einstein equations are quadratic in the dilaton gradients and form field strengths. 
Therefore, working to linearized order in the perturbations, the Einstein equation decouples 
from the perturbations of (/>, Hmnp, Fm and Fmnp since the latter do not have profiles 
in the AdS^ x background. Similarly the only contributions to the five-form equation 
of motion at linearised order are from the metric perturbations and the five-form field 
strength perturbations. 

The fluctuations can be expanded in harmonics [55]; 

hmn {x,y) = 

hma{x,y) = Y,iBll)Jx)YF{y) + Bll^Jx)DaY^^{y)) 
hiab){x,y) = + + 

ha{x,y) = ^TT^Hx)Y^^{y) (5.4) 

and 

fmnrst{^‘!y) — E I '■(!/) 

fmr,(x,y) = ^(bi;,nr,{x)OaY'' (y) + 

fabmnri^jy^ Y,mrnbl^r]i^)Yla\^iy) " 2bF^,{x)D^^Y^(y)) 
fabcmn{x,y) = 

fabcdm{x,y) = {x)eabc/DeY^^{y) + - 4)6^(x)e,,,/y/ny)) 

fabcde{Xjy) E bll){x)A^yabcdeY^^{y) (5.5) 

Numerical constants in these expressions are inserted so as to match with the conventions 
of [55]. Parentheses denote a symmetric traceless combination (i.e. = l/2{Aab + 

Aba) — l/Sfi'afe^a)- '^[ab] denote Scalar, vector and tensor harmonics 

whilst KF and A^® are the eigenvalues of the scalar and vector harmonics under (minus) 
the d’Alembertian. The subscripts t, v and s denote whether the field is associated with 
tensor, vector or scalar harmonics respectively, whilst the superscript of the harmonic label 
In derives from the number of components n of the harmonic. 

Not all fluctuations are independent - some are diffeomorphic to each other or to the 
background. This issue can be dealt with by imposing a gauge; for example, the de Donder- 
Lorentz gauge fixing condition is 

D'^h^ab) = D^ham = 0 (5.6) 

which sets to zero the coefficients ^(1)' ^ more elegant way of dealing with this 

issue is to construct gauge invariant combinations of the fluctuations. Such gauge invariant 
combinations of the fluctuations were constructed in [7], with the combinations reducing 
to the de Donder-Lorentz gauge fluctuations on imposing this gauge. 

In [7] the equations of motion satisfied by the fluctuations were constructed to quadratic or¬ 
der in the fluctuations, and the relation between five-dimensional fields and ten-dimensional 
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fields was also constructed up to quadratic order in the fluctuations. In the current context 
we are only interested in the five-dimensional Einstein metric and we work only to linear 
order in the fluctuations. We can therefore read off from [55, 7] the relationship between 
the five-dimensional Einstein metric perturbation and the ten-dimensional fields as 

uE _ 1^0 1^0 o ('c 

'"mn '^mn ^ ^ ^ ^mni V^* * / 

where the superscripts indicate that these are zero modes, i.e. associated with the trivial 
constant scalar harmonic. 

The type IIB supergravity equations lead to the linearized equation for the Einstein metric 

{Ce + = 0, (5.8) 


where Ce is the Einstein operator, defined as usual by 


^ + DpOmh^^ + DpDnh^ - DmDnhp^) . (5.9) 

The hve-dimensional equation of motion in turn follows from reducing the ten-dimensional 
action 


hiB = 


1 


ZKio 


—(lei{gMN) (r{9mn) — 


+ ... 


over the hve-sphere^. This results in 


d — 27r^ j ^ ^\/ det(^'mn) (^^R,{,9mn) + ' ‘ ‘ ^ , 


where we use the relation 


J-F _ _ 1 _ 

2K?n 2K?n 2 k| ’ 


(5.10) 


(5.11) 


(5.12) 


'"10 ^^10 

which is applicable when the AdS radius L is set to one. Thus the effective Newton constant 
is given by 

1 

Eor the probe brane system, the type IIB supergravity equations are solved with source 
terms, from the D-brane action, which in turn implies that the linearized Einstein equations 
in five dimensions are sourced. The complete ten-dimensional action is 


I = IlIB + Id7 


(5.14) 


where 

1^7 = -T7 J A J d^aS{x^ - X^{a^))e-^^-det{gMNdaX^dfiX^) + ■■■. (5.15) 

®Note that the lOd action must be supplemented with a self-duality constraint for the five form field 
strength. 
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Here a" denote the world-volume coordinates and the ellipses denote terms involving the 
world-volume gauge fields and Wess-Zumino couplings. The latter do not contribute in the 
case of the D7-brane embeddings under consideration here. 

The source term results in a stress energy tensor [56] 

tMN = ~ (5.I6) 

J V-det(5rMAr) 

where we denote the worldvolume induced metric as = guNdaX^dpX^ . The sourced 
IIB equation is thus 

{RmN - -^FmPQRsFn^^^ H-) = «lo(^MAr “ -^TgMN) = kIqTmN, (5-17) 

D o 

with T = g^^TMN- The trace adjusted stress energy tensor can be expanded in 
harmonics, using the same harmonic basis as for the metric: 

fmn{x,y) = Y,TL{x)Y^^{y) (5.18) 

%b){x,y) = + 

f:{x,y) = Y,f^^{x)Y^^{y) 

The correction to the five-dimensional Einstein equation only depends on the following zero 
modes (see Appendix A): 

{Ce + 4)/lmn = '^lo(^mn + 9mn) = ^mn, (5.19) 

where we have used the fact that the D7-brane embedding of interest does not source the 
RR five-form field strength. 

Given the D7-brane embedding, i.e. 9{z) = cos“^(m2:), the ten-dimensional energy mo¬ 
mentum tensor source can be computed as 

Tmn = -TrTMN6{0 - 0(^))<5(</>), (5.20) 

with 

Tzz = ^{1 - %u = (5.21) 

7^3 = (1 - m^z'^fgn^; 7^0 = 0 ; 

Tee = rn‘^z^{l — m^z^). 

where note that T = g^^TuN = 8(1 — vTz'^). Here we denote the metric on the unit 
three-sphere by go,^. The energy momentum tensor source can be projected onto spherical 
harmonics using Fourier decompositions of the delta functions: 



By projecting onto the zero mode, one can then immediately show that the ten-dimensional 
energy momentum tensor source is such that 

= r°, = 0, (5.23) 

with 

= -2), (5.24) 

TT^ 

and hence using (5.19) we find that 

izz =t^y = ^{l-rr?z^){m^z^-2)ri^y, (5.25) 

where as before 

to = IQtiGn{2tt‘^Ti). (5.26) 

Given the five-dimensional stress tensor tmn = imn — it is straightforward to see 

that the perturbation of the Einstein metric induced by this source is in agreement with 
that given in (3.14) and (3.18). 


5.2 Generalizations to other probe brane systems 

For a general probe brane system, the complete ten-dimensional action is 

I = IlIB + I Dp 

where 

I Dp = —Tp 


(5.27) 


:J J d^~^^a6{x^ — {a°'))e —det{'ycti3 + (5.28) 

+Tp j d^^x j 5{x^ -X^{o^)) ^^Cq 


with 


lap = OMNdaX^d^X^; (5.29) 

Xap = BMNdaX^dpX^ + F^0-, 

Ca,-a, = 

We consider embedding a brane into a type IIB background which is either AdS^ x or 
AdS^ Schwarzschild x 5^, so that the only background field profiles are for the metric and 
the five-form. Following the arguments in the previous section, we therefore only need to 
consider the equations for the metric and five-form perturbations, as the other perturbation 
equations decouple. 

The energy momentum tensor source is [56, 57] 

tMN = _T f ~ ^ ( 5 . 3 O) 
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where we define = 7 q /3 + TajS with being its inverse. The source in the five form 
equation of motion is 

Om = (5.31) 

2k\^Tp I dP+^a5{x^ - X^{a^))e^^-^^+^Fa,a, • • • 

Note that for a D7-brane this term only contributes if A 7^ 0. Therefore, provided 
that F A F = 0, the correction to the five-dimensional Einstein equation due to source 
D7-branes still depends only on the stress energy tensor zero modes: 

{Ce + = l^iQ{Tmn + = tmn- (5.32) 

One can thus compute the perturbation to the five-dimensional Einstein metric by project¬ 
ing the brane energy momentum source onto the appropriate combination of (spherical) 
zero modes. It would be straightforward to relax the condition F A F = ^ and obtain the 
correction to the five-dimensional Einstein equation, taking into account the sources in the 
RR field equations, but we will not analyse this case in detail here. 


The analysis above immediately allows us to treat D7-branes at finite mass, density (and 
temperature). The finite temperature background can be written as 

ds^ = —^dt^ + fdx^ -b -I-dfls; (5.33) 

= hp) = i + f,- 

with the temperature being T = '/2uqI'!i. The D7-brane embeddings can be expressed in 
terms of two scalar functions x{p) and a{p): 

9{p) = cos~^{x{p))] Fpt = dpa{p) = E{p), (5.34) 


where the potential is At = a{p). These embeddings can be found numerically, see [44, 45]. 
The main feature is that at any finite density, i.e. whenever the asymptotic form of the 


potential is 

4 d 

At — p - y -b • • • 


(5.35) 


with non-zero charge density d, the embeddings do not close off at finite radius. At finite 
temperature the embeddings have a spike which extends into the horizon, while at zero 
temperature this spike passes through the Poincare horizon. In other words, asymptotically 
as p —)• oo, the brane wraps the equator 9 = 7r/2 of the five sphere but there is a spike, 
0 —7- 0 as p —7 - uq. In the zero temperature limit, the spike solution becomes analytic for 
0 : 

E{p) ^ £■ 9{p) ^ 9ip, (5.36) 


with £ and 9i constant. 
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Focussing on the zero temperature limit for simplicity, the effective source stress energy 
tensor of (5.32) is given in terms of 6[p) and E{p) by 


tpp 

itt 

tij 


sin^ 6 ( 

-to -^-r - 

2p2(l + /,202_^2)2 

p‘^ sin^ 9 ( 

°2(1 + p202 _^2)| 3 

p^ sin^ 9 ( 

^°2(1+ p 2,92 _ ^2)i 3 


2 

3 





(5.37) 


Here 9 denotes dp9{p). To compare with the previous sections we change coordinates to 
z = \lp, and use the five-dimensional stress tensor tmn = imn — \tg^^: 


tzz - tf) 


sin^ 9 


2z‘^{l + z^9l - E ^)2 


ttt — to 


sin^ ( 


2(1 + ^202 -^2)5 


r(i + ^X^); 


tij — to 


2 • S 

z sm I 


2(1 + z202_^2)2 


j{l + z^9l - E"^) 5; 


ij, 


(5.38) 


where 9z = dz9. 

The metric perturbation induced by such sources can then be expressed as 

S{ds‘^) = + ^^t^^dx^dxi. (5.39) 

zZ ^z 

As previously the gauge invariant combination is 

f{z) = f{z) + zh'{z) (5.40) 

and there are now two independent Einstein equations: 

-/ - \ A9 - h') = ^zHzz] (5.41) 

^ 2 ;(/i' - g) + - h") = z^{tu + U), 

where we define tij = USij. These equations can be integrated to give 

Rz) = -^z'^tzz j J Rt + R ; 

dio 

dzf J + ti{w)). 

These equations can be solved analytically as z —)■ 0 and z —)■ oo. The near boundary 
expansions of the fields x and E are 


X = mz + cz^ -!-••• (5.43) 

E = 2dz^ H-, 
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where m is the quark mass, c determines the quark condensate and d is the density. The 
corresponding asymptotic expansions of the gauge invariant metric perturbations are 

f{z) = ^ + O mcz*) + •••); (5.44) 

{g{z) - h{z)) = *-^dz^ H- 

Note that in fixing a Fefferman-Graham gauge as z —)■ 0 one needs to take into account the 
shift in the AdS radius. The Fefferman-Graham gauge is obtained by choosing 

f{z) = (5.45) 

which then implies that 

zh'{z) = — {—nn?z^ + O (m‘^z‘^, mcz^) + •••), (5.46) 

6 

and hence 

h{z) ~ ^ “ w?z^ + O mcz*) +•••), (5-47) 

where the integration constant is fixed by the AdS radius. 

In the opposite limit of z —)• oo we can use the spike solution (5.36) to show that 


hz) 


{g{z) - h{z)) 


42 : 3(1 — 7 J 

42 ^ 3(1 _£: 2 )| ’ 


(5.48) 


and hence the metric perturbations are bounded in the deep interior. 


The effect of the metric perturbation (5.39) on the entanglement entropy is expressed in 
exactly the same way as in previous sections, since g{z) does not enter the entanglement 
entropy. Thus for a slab, following (3.22), the brane contribution to the entanglement 
entropy is 


6S 


L2 

4Gjv52 




fizs) 



(5.49) 


where z is the turning point of the original minimal surface and a = e/z. At zero density 
f(z) is zero for z > 1/m and continuity of the metric and its derivatives requires h{l/m) = 
h'{ljm) = 0. At any finite density /(z) is non-zero at finite z and there is no need to 
impose that the function h(z) vanishes at a finite value of z. In the finite density case the 
integration is therefore over the entire entangling surface, i.e. the upper limit 6 = 1 and 


SS 


AGnz'^ 




f{zs) + 





(5.50) 


since the other boundary term vanishes at s = 1. 
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Note that SS has no discontinuities in its derivatives with respect to mass or to the width of 
the slab at finite density since f{z) has no discontinuities in its derivatives at finite density. 
This provides another reason for viewing as unphysical the discontinuities discussed earlier. 

It is useful to define the difference between the entanglement entropy at finite density and 
that at zero density, for the same mass: 

r2 / ,/i _ o6 

= 4^ [i + 

where 

<5/(^) = f{z) - ^(1 - (5.52) 

Sh{z) = h{z) - h{z)^^^, 

where we have used the analytic expression for f(z) at zero density. From the asymptotic 
expansions (5.44), we can infer that the asymptotic expansion of 6f is 

Sf{z) = toO(mcz^) + 

We can also always choose a gauge such that 

dh(z) = toO(z^) + • • • ; (5.54) 

this simply corresponds to matching the gauge asymptotically at zero and finite density. 
Substituting into (5.51) the difference between the entanglement entropy at finite density 
and that at zero density is UV finite^ and only the integrated term contributes: 

<55 - <55j^o = ^ ds^^^^^<5/(zs). (5.55) 

Note however that this quantity does have discontinuities in its derivatives at z = 1/m, 
since the zero density quantity has such discontinuities. 


(5.53) 


6h{e 


Vl - a® 


(5.51) 


5.3 Numerical calculation of the entanglement entropy at finite density 


We now consider in detail the embeddings in (5.34) with the aim to explicitily carry out 
the computation of the flavor entanglement entropy in the case of finite density. We follow 
the analysis in [44, 45], though unlike the latter we focus on the zero temperature case. 
The background, instead of (5.33), is therefore AdS^ x 


ds^ 


p^[—dt^ + dxl] + 



+ dO^ + sin^6<<in| + cos^Od/)^ 


(5.56) 


and the probe D7-brane extends in {t, X 3 , p, Us}. We consider an embedding 9{p), and 
in addition we introduce a U{1) gauge field At{p) on the worldvolume of the D7-brane in 

^Earlier discussions of the UV finiteness of terms in the entanglement entropy induced by a chemical 
potential may be found in [58, 59, 60]. 


- 28 ~ 



order to study the gauge theory at finite density and chemical potential. The DBI action 
for this probe brane then evaluates to 

Id7 = -T7 j + P^{dpX? - 2(1 - X^)F% (5.57) 

where x{p) = cos [0[p)] and Fpt{p) = dpAt{p) is the electric field. The equation of motion 
for the gauge field has solutions with asymptotics given by (5.35), and since Id 7 does not 
depend explicitly on At^ there is a constant of motion d = 5lD7/dFpt. 


For solving the resulting equations of motion it is useful to eliminate the gauge field At 
from the action by performing a Legendre transform with respect to d. The equation of 
motion for y can then be obtained from the Legendre transformed action I^j as 


dp 


- X^)X L 

\/i - x' + V 


P^X 


1 + 






3(1 - + 2p^x^ 


(5.58) 


2A(F 


1 - + p‘^x^ 

/ 3®(1 — x'^)^ + 


where x = dpX and d = d/ToT- It is straightforward to show that asymptotically solutions 
to this equation take the form given in (5.43). 


We solve (5.58) numerically for a given d with regular boundary conditions imposed in the 
deep interior. Recalling that x = cos 6 the spike solution is such that x(0) —)■ 1. However, 
since xip) = —sin00, x(0) = 0 and is independent of the value of 0(0). Therefore we 
instead set boundary conditions at p = po 1: 

x{po) = 1 - x{po) = -a, (5.59) 

with a > 0 and (5^ <C 1. These boundary conditions correspond to 

0(po) = Hpo) = (5.60) 

Such conditions are consistent with the spike solution 0 = 0ip + • • • in (5.36) provided that 
a ~ S'^/po] if the latter condition is satisfied the solutions can be smoothly continued to 
p = 0. These boundary conditions differ from those used in [44, 45] due to the fact that we 
work at zero temperature. Note that the quark mass can be extracted from the embedding 
using limp^.oo(px) due to (5.43). 

The equations of motion for the gauge field are then given by Hamilton’s equations with 
as Hamiltonian, which reproduce the fact that d is a constant, together with the equation 

d,A = ( 5 . 61 ) 

This equation can be integrated to give 

AM = 2i j\p' , + ( 5 . 62 , 
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(a) J = 200, m = 200 (b) d = 0.5, m = 20 


Figure 5: Plots of (5/(z) for various values of d and m. In all cases d jiz) asymptotes to zero as 
z —>■ oo. 


where we have set At{p) —>■ 0 as /? —)• 0. The chemical potential p, is then given by ^ 44 ( 00 ) in 
the previous expression. Once the embedding x{p) has been found above, one can compute 
Fpt = dpAt and p from (5.61) and (5.62) respectively. Note that the parameter d indeed 
characterises the density since 

lim (p^dpAt) = 2d, (5.63) 

p^OO ^ 

in agreement with (5.35). 


To compute the entanglement entropy one must change coordinates as in the previous 
section. At zero temperature the coordinate transformation is trivial and simply amounts 
to setting z = l/p whilst leaving the other coordinates unchanged. One can then compute 
the five-dimensional stress tensor components in (5.38), and thus the gauge-invariant metric 
perturbation f{z) as defined by (5.42). We plot 5f{z), as defined in (5.52), in Figure 5 for 
various values of d and m, the latter being fixed by the choice of both d and x'(0). The 
same general features are observed for all values of the parameters; 6f{z) peaks around 
z = 1/m and has a long spike slowly asymptoting to zero as z —)> 00 . Although one might 
expect intuitively that the thickness of this spike is determined by the ratio d/m (with a 
larger ratio leading to a thicker spike), the results indicate that it is in fact the magnitude of 
X^(0) that determine this thickness, with a larger value of x^(0) corresponding to a thicker 
spike (and a larger magnitude of Sf{z) overall). 


It is a simple matter to compute the background subtracted entanglement entropy using 
(5.55): the result for d = 200, m = 200 is shown in Figure 6. The graph shows the 
entanglement entropy as a function of the depth of the entangling surface z, which is 
proportional to the slab width 1. It follows from (5.53) that the subtracted entanglement 
entropy increases quadratically with z for z <C 1; if d/ = Az^ then 


6S 


. L^X v/^r(I/3) ,2 

48Gn r(Il/6) ^ 


(5.64) 
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Figure 6: Plot of the background subtracted entanglement entropy for d = 200, m = 200 as the 
width of the slab is increased. The entanglement entropy increases quadratically until the depth 
of the entangling surface is 1 jm and then slowly saturates to a constant value as the width of the 
slab is increased further. 


The metric perturbation bf reaches a maximum around z ~ 1/m and is very small for 
z > 1/m, and therefore the entanglement entropy of surfaces which extend to turning 
points z 3> 1/m saturates. 

The D3-D7 system has a rich structure of phase transitions as the chemical potential, 
temperature and magnetic field are varied, see [61, 62]. It would be interesting to use 
entanglement entropy to explore these phase transitions, extending the above results. Note 
that the entanglement entropy for massless flavors at finite density was discussed in [12]; 
our method would give the same results for massless flavors, and it would be interesting 
to explore how the entanglement entropy changes as one increases the ratio of density to 
mass. 


6. Field theory interpretation 

The D3-D7 system is dual to AA = 4 SYM coupled to AA = 2 massive hypermultiplets. 
The key features of this field theory are as follows. The field content of AA = 4 consists 
of the gauge fields scalars transforming in the fundamental of the R symmetry 
group SO{6) and spinors A® transforming in the spinor representation of S'0(6). The 
hypermultiplets consist of scalars y and fermions r] transforming in the bifundamental 
of the SU{Nc) and SU{Nf) gauge groups. In the massless case the addition of these 
hypermultiplets preserves an 50(4) x 50(2) subgroup of the R symmetry group of A/" = 4 
SYM. The hypermultiplets are coupled to the AA = 4 SYM fields by potential terms of the 
form 

1 = 1 d\TTsu(N.)TTsuiNf) ( 6 . 1 ) 

Separating the branes by a distance m (in string units) corresponds to introducing a mass 
term m for the hypermultiplets, which breaks the conformal invariance and breaks the R 
symmetry group further to 50(4). 
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There are two distinct regimes of interest: the mass parameter m being small relative to 
energy scales of interest, and the mass parameter being large compared to scales of interest. 
In the former case the theory is clearly described in terms of a small mass perturbation 
of a conformal field theory, and we can use the underlying conformal invariance to un¬ 
derstand the entanglement entropy. Entanglement entropy for relevant perturbations has 
been studied recently [23] and we will discuss the relation to our results below. 

In the opposite regime, of high mass, at energy scales much lower than m we can integrate 
out the hypermultiplets, effectively setting X ~ The potential term above controls the 
leading deformation to the N = A SYM theory: at low energies the effective description 
must be 1 /■ 

I = IsYM H-2 / (6.2) 

J 

where Oq is an operator of dimension six in the J\f = 4 SYM theory, i.e. it is an irrelevant 
deformation of the SYM conformal field theory. The dimension six operator explicitly 
breaks the R symmetry from to 50(4) and is therefore charged with respect to the 

50(6) R symmetry of AA = 4 SYM. The backreaction of this deformation on the stress 
energy tensor, which is an R symmetry singlet, is necessarily quadratic in this deformation, 
i.e. the stress energy tensor is only affected at order 1/m^. The behaviour of entanglement 
entropy under irrelevant deformations has been less studied and we will explore this case 
in more detail below. 


6.1 Zero mass: marginal deformation of CFT 

In the limit of zero mass, the brane contribution to the entanglement entropy of the slab is 

tpL^ ( 3 , 0f r(-l/3) ^ 

4 SGn V2e2 ^ 2^2 r(l/ 6 ) J ’ 

where implicitly we have fixed a gauge choice such that 


h(z) = —. 
^ ^ 12 


(6.4) 


Note that this expression is proportional to the entanglement entropy of a slab in AdS^: 

(6.5) 


tpL^ ( 1 V^ T{-l/3) \ 

2 Gn 6^2 r{l/6) J ' 

The entanglement entropy of a spherical surface in AdS^ is 


TT 




1 




( 6 . 6 ) 


The brane contribution to the entanglement entropy for a spherical surface at zero mass is 


6S = 


tpir / i?2 1 


8Gn V2e^ ' 2 
which is again proportional to the AdS result 


+ -log{€/€lR) - 


1 


(6.7) 
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This is to be expected: suppose that the AdS radius is scaled as 

LAdS ^AdS (1 + ^Ads) ( 6 . 8 ) 

with 6AdS ^ 1. Since the bulk entangling surface is of dimension three, this implies that 

S^S{l + 5Adsf^-s(^l + ^SAdsy (6.9) 

In the case at hand, the effect of the brane is to shift the AdS radius as 

^Ads = (6.10) 

The brane contribution to the entanglement entropy at zero mass is therefore precisely 

-55 = l^AdsS, ( 6 . 11 ) 

explaining the results above. 


6.2 Small mass: relevant deformation of CFT 


In the small mass regime, the mass m is smaller than the energy scales of interest, so the 
system can be viewed as a mass perturbation of an underlying CFT. 

Half space: Let us consider first the case in which the size of the slab Z —)■ oo, i.e. the 
space is divided into two regions by a plane; the brane contribution to the entanglement 
entropy is given in (3.31). In particular the logarithmic divergence is 

<55 = yy^ m^log{me) = ^Tq{w?L^) log(me), (6.12) 

where in the latter expression we use (3.19) i.e. Tq is the effective tension of the D7-brane, 
reduced over the three sphere. 

The mass deformation in the field theory is associated with the following slipping mode of 
the D7-brane on the three sphere: letting the five sphere metric be 

= dO"^ + sin^ + cos^ Odcj)^ (6.13) 


then the slipping mode is associated with the angle 6, i.e. we retain only the following 
terms in the D7-brane action 


I = To 


d^x^/gsin^ 9^/1 + 


(6.14) 


where we have integrated over the three-sphere, because the mode of interest is an 50(4) 
singlet and hence there is no dependence on the three sphere coordinates. The metric g^y 
denotes the AdS^ metric and we work here in Euclidean signature as we need to compute 
correlation functions. The resulting equation of motion for 9 is 


0 = D^-dcot^ 


1 g^^^d^9{gp-dp9d^9) 

2 1 -h g^^''^^9^u9 


(6.15) 
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where □ is the Laplacian in the Euclidean AdS^ metric. Linearising this equation around 
9 = 7^12 gives 

0 = 061+ 36*, (6.16) 

i.e. the scalar is dual to an operator of dimension three. 

As we will discuss shortly, we are interested in computing the normalization of the two 
point function of this operator, and it thus suffices to consider the solution to the linearized 
equation of motion (6.16). From [27], we can read off the operator one point function in 
terms of the asymptotic expansion of the scalar field: 

(Oa) = To ^—20(2) + -0^0) + '^(0)^(0)^ (6-17) 

where 

0 = 2;(0(o) + ^(1)2^ + ^(2)2^ + ^(2)2^ log(2:) + • • • ) (6.18) 

Here □(q) refers to the Laplacian in the boundary metric ^((o) and ii(o) is the scalar curvature 
of this metric, which is zero in our case. Working in momentum space, the regular solution 
to the linearized field equation (6.16) is 

e = ei^){k){z^Ki{kz)), (6.19) 

with k the momentum and 9(^o'j{k) corresponding to the Fourier transform of the source. 
One then expands (6.19) about z = 0 to identify the various terms in (6.18), for example 

^{2){k) = 6'(o)(/c) Q(-l + 27 )A :2 + ^fc^log (6.20) 

where 7 is the Euler constant. Functionally differentiating the one point function (6.17) 
with respect to the source 0 ( 0 ), and setting the source to zero, therefore gives 

{0^{k)0^{-k)) = ne log (^) + • • • > (6-21) 

where contact terms have been dropped. Fourier transforming back to position space gives 

(O3(x)O3(0)) = , (6.22) 

where TZ denotes differential renormalization, see [63]. 

In [23] it was shown that for a CFT deformed by a relevant operator 


I + A / 


(6.23) 


the change in the entanglement entropy of the half space is 


/ 7 \ d~|~2 y V 

SS = Lb ^log (55:) , 

4(d-l)r('^) 


(6.24) 
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where ejjv and em correspond to UV and IR cutoffs, respectively, A is the area of the 
dividing surface while N is the normalisation of the two point function of O, i.e. at 
separated points 

{0{x)0m = (6.25) 

with A = {d + 2)/2 the (relevant) operator dimension. Note that the entanglement entropy 
is unchanged to first order in the perturbation A. 

In our case the normalisation is given in (6.22), X = m and A = L^. Hence 

iS = log (, ( 6 , 26 ) 

which exactly agrees with (6.12), taking the IR cutoff to be ejn = Ijm. 

Sphere: For the spherical entangling region the logarithmically divergent terms are 

^ 0 log(me) = To log(me). (6.27) 

The second of these terms was explained above. The first is proportional to the mass 
deformation and can be expressed in the same form as (6.24), with 

A'T' 

A = iTrR^-, M=^, (6.28) 

setting the IR cutoff to be = 1/m. 

Slab: For a slab of finite width, the logarithmically divergent terms are precisely twice 
those for the half space: 

SS=?^m^L^kg(‘XL], (6.29) 

Taking into account that the area of the entangling surface is in this case 2L^, we again 
find exact agreement with (6.24). 

At first sight it may seem surprising that the expression (6.24), which was derived for 
the half space using the known modular Hamiltonian, is applicable to other entangling 
geometries, with the entangling area replaced by the appropriate value. However, from the 
field theory perspective, one could derive the result for the spherical region by conformal 
transformations (the modular Hamiltonian is also known) and the divergent contributions 
for a slab, being local, must necessarily give exactly twice the result for an infinite slab. 
It was argued in [23] that the result should hold for any geometry, since any entangling 
surface is locally flat; of course for a curved surface there are additional contributions 
to the entanglement entropy beyond this universal contribution. The dilaton effective 
action approach was also used in [64, 65] to derive the logarithmic divergences for any 
shape entangling region, up to a universal coefficient computable from the dilaton effective 
action. 
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There is also a very simple holographic way to understand why the formula (6.24) is 
applicable to the logarithmic divergences of any shape entangling surface, generalising the 
work of [66, 26]. Deforming the conformal field theory by an relevant scalar operator 
corresponds in the bulk to coupling gravity to a massive scalar <1> i.e. we consider 

I = J (^R + did + 1) - 2 + •••), (6.30) 


where = A(A — d), with A < d the dimension of the dual operator. Here implicitly 
we are working perturbatively in the scalar field so we include only quadratic terms in 4> 
with the ellipses denoting higher order terms. The normalisation of the operator two point 
function is [67, 63] 


(O(x)O(O)) 


( 2 A - d)r(A) 

16 TrGi\[TT 2 T{d — 


for A = d /2 + k with k an integer. In particular we can 


1 

A'l ’ 

2 1 

write 


(6.31) 


r{d/2 +1) 

g^d/2+lQ 


for A = d /2 + 1. 

Working perturbatively around an AdSd+i background, a scalar field profile 


$ = 


(6.33) 


corresponds to deforming the field theory by the dimension A operator, with A charac¬ 
terising the deformation. At quadratic order in the source there is a backreaction on the 
metric. Letting the metric perturbation be as before 

5ids‘^) = ^ [fiz)dz‘^ + hiz)dx^dXfj/) (6.34) 

then the Einstein equation implies that the gauge invariant combination of these pertur¬ 
bations is given by 

fiz) + zh'iz) = (6.35) 

Working in Fefferman-Graham gauge we may set /(z) = 0 in which case 

hiz) = - ^ (6.36) 

4(a — i) 

Now consider an entangling surface in the deformed metric. Let the induced metric for the 
minimal surface be and fix a static gauge such that 

Z = z; = cj“; A* = X\z, a“) (6.37) 

where asymptotically as z —)• 0 

A*(z,cj“) = A*(a“) + --- (6.38) 
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The entanglement entropy contains divergent terms from 

= 4 ^ / "" 4 ^ / + • • • + + ■■■) ( 6 . 39 ) 

where 

rab = daX\a'^)dbX^{a'^)S,j ( 6 . 40 ) 


is the induced metric for the entangling surface on the boundary. Integrating over the 
radial coordinate one finds the usual power law volume divergence 

4G^,(<i - 2 ),“-^ / . (6.41) 

with z = e being the UV cutoff and A being the volume of the entangling surface in the 
boundary. From the relevant perturbation one obtains a logarithmic divergence whenever 

A = i(d + 2) (6.42) 

where 

(55 =-—(d-2)A^/iolog(e/e/ij). (6.43) 

OVrAT 

Using the identity (6.32) we thence obtain 


(55 


AfAX'^ 


{d - 2)7r‘^/2+i 

4(d- l)r(d/2 + 1) 


log(e/e/R), 


(6.44) 


which is exactly (6.24) but does not assume any geometry for the entangling surface. It 
would be interesting to prove this result from the field theory; the modular Hamiltonian 
for generic entangling surfaces is not known but the holographic result suggests that one 
should be able to compute the logarithmic divergences without complete knowledge of the 
modular Hamiltonian. 


Note that this result has a straightforward generalisation to irrelevant deformations. De¬ 
forming the conformal field theory by an irrelevant scalar operator corresponds in the bulk 
to coupling gravity to a massive scalar with = A(A — d), with A > d the dimen¬ 
sion of the dual operator. Following the same steps we see that the entanglement entropy 
contains divergent terms from 

® = 4^ / 4^ / pT(l + + ...) (6.45) 


but since A > d the metric perturbation h{z) always gives rise to additional UV divergences: 


(55 = 


A2(d-2) ^ 

32GAr(d - 1)(2A - d - 2) ' 


(6.46) 


Using the identity (6.32) we thence obtain 


4^x2 (d-2)7r^/2+l _^ 

4(d - 1)(2A - d - 2)F(d/2 + 1) e2A-(i-2 ’ 


(6.47) 


where Af is the operator normalisation. As for the usual power law divergences, such terms 
are not universal but nonetheless will be given an interpretation in the following section. 
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6.3 Large mass: irrelevant deformation of CFT 


By large mass, we mean that the mass scale m is higher than the energy scales of interest. 
This implies in particular that m ^ 1/1, where I characterises the size of the entangling 
region, i.e. the width of the slab or the radius of the sphere. Thus we are always working 
in the regime fi ^ 1. The UV divergent contributions to the entanglement entropy have 
already been explained above and here we are interested in explaining the leading finite 
contributions for ^ ^ 1. To decouple such contributions from the divergent terms, it is 
useful to look at the differentiated quantities (3.51) and (4.23) (which we argued previously 
do not receive contributions from UV divergent terms). 

As stated above, for large mass, the effective IR description is in terms of an irrelevant 
deformation of SYM. It is easy to understand the effects of such a deformation on the 
entanglement entropy from the dual perspective. For a deformation by an operator of 
dimension six with A ~ 1/m^ the change in the metric behaves as l/{mz)'^. The metric 
perturbation can only be viewed as small relative to the background AdS^ metric when 
(mz)^ » 1. 

Slab: Now it is straightforward to infer the effect of the irrelevant deformation on the 
entanglement entropy of a slab. Since the latter scales extensively with the volume of the 
slab, L^, and the metric is corrected at order 1/m^, the change in the entanglement entropy 
goes as I?/m^. The underlying theory is conformal so the only scale in the problem is 
the width of the slab 1. Since the entanglement entropy is dimensionless the effect of the 
massive modes is to change the entanglement entropy as 


55 ~ 




(6.48) 


which indeed agrees with the term found in (3.51). 

We can infer this answer using (6.47), which in the case of A = 6 and d = 4 gives 


55 ~ 



(6.49) 


In the case at hand A ~ l/nr?. The description in terms of an irrelevant deformation of 
SYM is only valid provided that we consider entangling surfaces for which I ^ 1/m. The 
effective cutoff should therefore be e ~ f ^ 1/m and hence we reproduce the formula above. 

Sphere: For a spherical entangling surface the leading contribution to the differentiated 
entanglement entropy at large ^ behaves as (4.23) 


55. 


m?R^ 


+ 0 


1 




(6.50) 


Thus although there is a 1/m^ term (as above) this is not the leading contribution. A 
simple way to understand the origin of this term is by exploiting the CHM map [2]. The 
entanglement entropy for the spherical region is then computed by computing the entropy 
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in the mass deformed theory on a hyperbolic space. Since the fields are conformally coupled 
the action contains the terms 

/ = ^ y {{dxf + > (6-51) 

where x is a hypermultiplet scalar and 7^ ~ 1 /E? is the Ricci scalar. Now when we integrate 
out the hypermultiplets we obtain additional terms in (6.2): setting x ~ 1/m we obtain a 
contribution from the curvature coupling of order which is in agreement with 

the expression above. 

6.4 Conformal perturbation theory at higher orders 

The brane contributions to the entanglement entropy for the half space (3.31) consist of only 
the divergent term arising from the shift in the AdS radius and the logarithmic divergence 
discussed above. Since the logarithmic divergence is expressed entirely in terms of the 
coefficient of the two point function for the dimension three (fermion mass) operator, this 
contribution to the entanglement entropy is trivially not renormalized relative to the weak 
coupling result: since the operator dimension is protected, there is such a contribution 
regardless of the coupling. The result (3.31) only includes powers of m up to m?, which 
follows on dimensional grounds: the dimensionless entanglement entropy scales extensively 
with the slab area and therefore the only way that contributions at order or higher 
could arise in the entanglement entropy would be if the latter was IR divergent, i.e. the 
contributions would have to scale as 

5S ~ (6.52) 

where the IR cutoff A/ij ;§> 1 (relevant perturbations cannot introduce UV power law 
divergences so an IR cutoff is the only possibility). However, the entanglement entropy is 
an IR safe quantity and therefore no such dependence on an IR cutoff should arise. 

From the perspective of conformal perturbation theory, it is not obvious that there are not 
contributions to the entanglement entropy from higher order terms, i.e. terms of order 
or higher. The change in the entanglement entropy is in general expressed as [21, 22, 23] 

5S = -m{OK) + {{OKK) - (OO)) + 0{m^) (6.53) 

where O is the deforming operator and K is the modular Hamiltonian. The first term 
vanishes by conformal invariance and the second term gives the logarithmic divergence. By 
the argument above, higher order terms (dependent on higher order correlation functions) 
must vanish and it would be interesting to show this explicitly. 

7. Differential entropy 

Given the expression for the entanglement entropy of a slab, we now proceed to compute 
the differential entropy. Following [31, 32, 33, 34, 35, 36], the differential entropy should 
correspond to the area of a surface in the backreacted geometry - see Figure 7. 
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z = 0 


Figure 7: An illustration of the equivalence between the differential entropy of a boundary partition 
and the area of a corresponding hole in the bulk. As the number of strips partitioning the boundary 
tends to infinity, the turning points of the associated minimal surfaces (blue) form a smooth hole 
in the bulk (red) whose area equals the differential entropy. 


The differential entropy is defined as 

OO 

E = J2[S{Ik)-S{hnIk+i)] (7.1) 

k=l 


where {Ik} is a set of intervals that partitions the boundary. We will cover the boundary 
with n intersecting slabs - we take Ik to be a slab of width Ax, and the intersection IkCiIk+i 
is thus a strip of width Ax — Lxjn where Lx is the regularised length of the x-direction. 
At the end we will take the limit n —>■ oo. 

For a slab in AdS^ there is a relation between the strip width Ax, and the maximum bulk 
depth of the associated extremal surface z, see (3.12). When the slab lies in the perturbed 
geometry this relation is modified (see for example (3.36)), and therefore it is useful to 
leave the relation implicit as 


Ax = cz = (co + Cl to + 0{tl))z. 


(7.2) 


with 



(7.3) 


and Cl depends on whether /r is greater than or less than one. The differential entropy 
takes the following form 


E = lim n[S'(Ax) — S{Ax — Lx/n)] 

n^oo 


(7.4) 
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where the overall factor of n arises from the fact that all slabs are of equal width. For 
AdS^ this gives 


E = 


-^ lim n 

12r(g)GAr 



(Ax - L^jnf 


(7.5) 


and hence 


6r(i)Gjv (Ax)3’ 


(7.6) 


where V = 
rewritten as 


L^Lx is the regularised three-volume. Using (7.3) this expression can be 


E = 


_^ 

AGn (Ax) 3 ’ 


(7.7) 


and the latter is manifestly equal to the volume of the turning point surface (i.e. the hole) 
divided by AGjy. 


7.1 Large mass 1 

Proceeding to calculate the flavor contributions, we begin with the very large mass case. 
The relevant expression for the entanglement entropy is (3.30). The two leading terms are 
manifestly independent of z and thus of Ax which implies that the differential entropy 
vanishes in this limit: 


5E = lim n[(55(Ax) — 6S(Ax — Lx/n)] —)■ 0. 


(7.8) 


The limit /x ;§> 1 means z ^ 1/m and thus the hole formed in the bulk is far away from 
the probe brane (which ends at z = 1/m). In the large jj, limit the leading contribution to 
the differential entropy arises from the third term in (3.30) and hence 


SE 


2304m4Gjvn-^” 

tpVct 

38Am^{AxyG]\f ’ 



Ax® 


(Ax - Lx/ny\ ’ 


(7.9) 


where cq is defined by (7.3). Note that since the gauge dependent terms are independent 
of Ax they also automatically cancel in the differential entropy. 

Now let us compare the differential entropy with the volume of a hole of radius z. Since 
the backreaction on the AdS^ metric in this region is zero one might naively expect that 
the change in the differential entropy is zero. However, this does not take into account the 
fact that the relation between the turning point radius z and the width of the boundary 
slabs Ax is modified (3.36). In other words, the gravitational entropy of the hole remains 


E 

4GAr z3’ 


(7.10) 
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since the metric is AdS^, but to express this quantity in terms of Ax we need to use the 
relation in (7.2) to first order in the perturbation, resulting in 


, ^ V f eg 3cgcifo \ ^ Vcl toV4 

gray \^(Ax)3 (Ax)^ / 4GAr(Ax)3 384m^(Ax)’^GAr ’ 


(7.11) 


where we use 


288m^5^ 288m^(Ax)^ 


(7.12) 


where the latter equality is to order to. The second term in (7.11) is in exact agreement 
with the calculation of the perturbation in the differential entropy (7.9); we will see in 
section 7.3 that this agreement holds for all ^ > 1. 


7.2 Small mass n < 1 


We proceed with the small mass case. Extracting from (3.28) the terms which depend on 
z one obtains 


5S 


( \/^r(-i/3) 1 r(l/3) 4 2 

48Gjv V 12 r(7/6) ^ 12 r(ll/6) 


2irr\ogz + 


(7.13) 


From these terms in the entanglement entropy we can now proceed to compute the differ¬ 
ential entropy using (7.4). Again it is clear that all divergent and gauge dependent terms 
cancel from the differential entropy, since they are independent of Ax. Noting that: 


(Ax)^ — (Ax — Lxlnf = 2Ax— -|- C>(l/n^) 

n 


1 


1 


2 L, 


(Ax)2 (Ax — Lx/nY Ax^ n 
log 


+ 0{1/tY 


/^x — Lx/n 1 Lx , 9 , 

- ^ = -w-^ + 0{l n^) 

Ax Ax n 


we obtain the perturbation in the differential entropy; 


toV ( 3cq 4 \/^ r(l/3) Ax 2m?\ 
48Gjv y ~ 2Ax3 “ ^ T r(ll/6)^ ^ ^J’ 


(7.14) 

(7.15) 

(7.16) 


(7.17) 


where cq is as given in (7.3). 

Now let us compare this expression with the gravitational entropy of a hole of radius zt- 

where we take into account the metric perturbation relative to AdS^. In this case the 
relation between zt and A is given by 

Ax = coZt{l - ^h{zt)) + c'ltozt, (7.19) 
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where cq is defined in (7.3) and 


/ 

Cl = - 


1 



Co + 




3 


9 r(f) J 


(7.20) 


Thus we note that the gauge dependent quantity h{zt) cancels from the entropy of the 
hole, with 

^grav = 4^^ (Ax) 3 ^'^oCl^o) j (’7-21) 

and moreover the change in gravitational entropy of the hole agrees exactly with the change 
in the differential entropy (7.17). 


7.3 // > 1 


We now proceed with the case ^ > 1, without requiring ^ ^ 1. The entanglement entropy 
for ^ > 1 is given by equations (3.29) and (3.30); 

(i + C 1). {2. 2), l/A*') - Y^Fi (-1/2, -1/3,2/3,1/,.®) 

+ ^2^1 (-1/2,1/3,4/3,1/^®) + 2fiHog{fia)] + SSscheme{m,e) 

(7.22) 


Recalling that fi = mz and a = ejz, the divergent and log parts are independent of 5 (and 
thus Ax) and so do not contribute to the differential entropy, nor do the gauge dependent 
terms. Also recall that z = Ax/c, where throughout the following we can replace c —)• cq 
to order Iq since 5S already contains an overall factor of to- We thus have: 


5S = 


48G VOm^Ax® 


({1/2,1,1). {2,2}. ^^) - (-1/2, -1/3,2/3, 


c® 

m^Ax® 


m 


+— 2 F 1 -1/2,1/3,4/3, 


m^Ax® 


+ ... 


(7.23) 


where the ellipses denote ^-independent terms. We want to compute the differential en¬ 
tropy: 


5E = lim n[(55(Ax) — 6S{Ax — Lx/n)] 

n^oo 


(7.24) 


which requires the relations: 


(-1/2,1/3,4/3, /,) - .F, (-1/2, l/3,4/3, 
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(7.26) 
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-. 3^2 { 1 / 2 , 1 , 1 }, { 2 , 2 },^ - 


X®/ {x — Lx/nY^ 


3F2 {1/2,1,1},{2,2}, 


(x - La,/n)6 


12L^ 


ax 


-1 + Wl - ^ ) - + 0{l/v?‘) 


X® / n 


(7.27) 


where here a/x® < 1 since we are considering /x > 1. For the former two results we have 
used the standard identity: 


-^ 2 Fi(a,b,c,x) = —2^1(0 + 1,6+ l,c + l,x) 
dz c 


and for the latter result we have used the following relation: 


3F2({1/2,1,1},{2,2 },x) = - log 

X 


1 + v/r^ 


— \/l — X + 1 


(7.28) 


(7.29) 


where, again, here x < 1 since we are considering ^ > 1. The differential entropy is then 
calculated to be: 


6E = 


toV q 


384G (Ax)3^4 
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12 5 1 


16^ -1 + A/l - -fi + 62T1 7:, 77,77, -F + 32^1 -, -, -, ^ 
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14 7 1 


2 ’ 3 ’ 3 ’ /x® 
(7.30) 


where recall V = L^L^, and we leave some Ax-dependence implicit in ^ for notational 
clarity. 

We now want to compute the corresponding change in the gravitational entropy of the 
hole. As shown in section 7.1, the change in gravitational entropy of the hole is given by: 


rp _ ^ Scgcito 

()^grav-4^^ Ax3 

where ci is defined by equation (7.2) and is given by equations (3.34)-(3.35): 
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(7.31) 


Cl = 
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12 5 1 


16/X -1 + W1 - ^ + 62 T 1 77> F’ F’ TT + Tx 
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14 7 1 


2 ’ 3 ’ 3 >® 


(7.32) 

The equality of the differential entropy and the gravitational entropy of the hole is then 
manifest. 


8. Entanglement and differential entropy for top-down solntions 

The main focus of this work has been to develop a systematic method for computing the en¬ 
tanglement entropy for top-down brane probe systems. Our method however immediately 
generalises to any top-down solution which can be viewed as a perturbation of AdS^ x 5® 
(or indeed a perturbation of the finite temperature AdS^ Schwarzschild x 5®): using the 
Kaluza-Klein holography dictionary we can extract the effective five-dimensional Einstein 
metric and thence compute the entanglement entropy. 
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It is worth noting that given a general asymptotically AdS^ x supergravity solution 
one cannot extract the five-dimensional Einstein metric, except asymptotically around the 
conformal boundary. Brane probe systems in which the perturbations relative to AdS^ x 
are small everywhere, and thus the Kaluza-Klein dictionary can be used to compute the 
five-dimensional Einstein metric at all scales, are special cases. In general the Kaluza-Klein 
dictionary becomes intractable when the metric perturbations relative to the AdS^ x 
background are of order one, i.e. at some finite distance from the conformal boundary. 

Another subset of ten-dimensional supergravity solutions can be expressed in terms of 
the uplifts of five-dimensional gauged supergravity solutions, for which the Einstein met¬ 
ric is known from the lower-dimensional theory. For example, specific cases of Coulomb 
branch solutions are realised as solutions of AA = 8 gauged supergravity in five dimensions; 
generic Coulomb branch solutions are however not realised as lower-dimensional solutions. 
The Coulomb branch examples also illustrate the fact that the causal structure of the 
five-dimensional Einstein metric is generically not the same as that of the uplifted ten¬ 
dimensional metric: the former can have naked timelike singularities which correspond to 
harmless null horizons in the uplifted solutions. One could thus envisage a scenario where 
the lower-dimensional metric had no entanglement shadows of the type discussed in [34] 
but the uplifted solution had shadow regions which could not be probed by five-dimensional 
fields at all. 

The entanglement entropy is computed from the five-dimensional Einstein metric and the 
differential entropy (built from entanglement entropy) therefore reconstructs areas of holes 
in the five-dimensional Einstein metric. The ten-dimensional metric cannot be recon¬ 
structed just from the five-dimensional Einstein metric: the uplift map requires all the 
matter fields in the lower-dimensional theory. Therefore the standard entanglement en¬ 
tropy and differential entropy cannot in principle reconstruct the ten-dimensional geometry 
without additional information. 


8.1 Coulomb branch examples 


To illustrate the above discussions we consider a particular Coulomb branch solution dis¬ 
cussed in [68, 69]. The Coulomb branch of AA = 4 SYM corresponds to the spontaneous 
breaking of the gauge symmetry by giving VEVs to the scalars - on the gravity side, these 
solutions are represented by multi-centre D3-brane solutions. These flows break supercon- 
formal invariance but preserve sixteen supercharges. 

[69] studies particular Coulomb branch solutions which admit consistent truncations. These 
flows are described in five-dimensional gauged supergravity by a single scalar field x(r) 
where r is a radial coordinate for the 5-dimensional metric 

ds^ = + dr^ ( 8 . 1 ) 


where r —)• oo corresponds to the conformal boundary. The BPS conditions can be written 
as 


^ ^ gc^ dA ^ 

dr 2 dx dr 3 


( 8 . 2 ) 


-45 - 





Figure 8: Plots of po and AS as functions of ^ for cr = 0.1. Both functions asymptote to zero as 
I oo. 


where W(x) is the superpotential and g is the gauged supergravity coupling constant. In 
particular, we will be interested in the solution which, from the ten-dimensional point of 
view, corresponds to the D3-branes being uniformly distributed on a disc of radius a in 
the transverse space, preserving S'0(4) x SO{2) of the SO{6) symmetry in the conformal 
AdS^ X solution. In this case one has 


W{x) = —— -e 


^(x) = ^log 




+ log (<t) 


1 _ eV6x 

where we have as usual set the radius of curvature of AdS^ to be one. 
Redefining the radial coordinate one can write the metric as 


ds^ = 


dx,dx^ + 


- ( 1 + ^ 


(8.3) 

(8.4) 


(8.5) 


where p —)• oo at the conformal boundary and the metric is AdS^ for a = 0. The uplifted 
ten-dimensional metric is then expressed as 


ds^Q = A + dsj^, 


( 8 . 6 ) 


where the warp factor A depends on the sphere coordinates and dsj^ is a metric on a 
warped sphere. Explicitly 

^ C = (1 + ^cos^e); (8.7) 

dsj^ = ^ (C^d0^ -|- cos^ 9dQ^ + A® sin^ Odcj)^) . 

The five-dimensional metric has a naked timelike singularity at p = 0 but the uplifted 
geometry has a null horizon at p = 0. 
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The five-dimensional metric satisfies an a-theorem; the warp factor A{r) in (8.1) decreases 
monotonically as r decreases. Correspondingly the entanglement entropy and the differ¬ 
ential entropy monotonically decrease as the scale of the entangling region is increased. 
Working for convenience in the coordinate system (8.5) the entanglement entropy of a slab 
is 

Q - W C + 

where A is the UV cutoff and po is the turning point. The quantity 

k2j2 

AS = S-— (8.9) 


is UV finite by construction. The relation between the width of the slab I and the turning 
point is 


poo 

I = + / 

j po 


dp 

cj2)(p4(p2 + ^2) _ pi(^p2 ^2)) 


( 8 . 10 ) 


The quantity AS{1) monotonically decreases as I is increased. We can express the differ¬ 
ential entropy associated with strips of width I as 


E{1) = 


d{AS) 

dl 


( 8 . 11 ) 


and this quantity also decreases monotonically with 1. This had to be true since (by 
construction) E{1) can be expressed in terms of the warp factor in (8.1), which is known 
to satisfy the a-theorem, see [68]. It is interesting to note that the differential entropy 
is proportional to the finite entropy (3.49), which in turn is known to play the role of a 
c-function in two spacetime dimensions. In any holographic background dual to a RG flow, 
the differential entropy is by construction expressed in terms of the warp factor in (8.1), 
which satisfies the a-theorem provided that appropriate energy conditions are imposed on 
the bulk stress energy tensor [68], and therefore the differential entropy has the correct 
property to correspond to an a-function (in any dimension). 


Entanglement for Coulomb branch geometries has been discussed in earlier papers [70, 71] 
from a ten-dimensional perspective, i.e. minimal surfaces in the ten dimensional geometry 
were explored. As emphasised throughout this paper, the standard entanglement entropy 
can only be computed from the five-dimensional Einstein metric, which can only be ex¬ 
tracted near the boundary in the case of separated brane stacks discussed in [70, 71], since 
no consistent truncation to five-dimensional supergravity exists (Kaluza-Klein holography 
allows us to extract the Einstein metric near the conformal boundary, where the metric 
is close to AdS^ but deep in the interior the metric is not close to AdS^ and therefore 
the method cannot be used there, see [7, 72]). Codimension two minimal surfaces in the 
uplifted geometry (8.6) compute the generalized entanglement entropy defined in [70]. It 
would be interesting to understanding the field theory interpretation of this quantity better; 
see recent discussions in [73]. 
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9. Conclusions 


One of the main results of this paper is a systematic method to compute the effective 
lower-dimensional Einstein metric for top-down brane probe systems, using which one can 
extract the entanglement entropy. We have illustrated this method with the case of quarks 
at finite mass and density, at zero temperature. It would be interesting to explore the 
finite temperature behaviour of the entanglement entropy and how it captures the phase 
transitions found in [44, 45]. Our method is applicable to any brane probe system and 
could for example be used to evaluate the entanglement entropy in models of quantum Hall 
physics [74, 75] and in top-down models of the Hondo effect [76]. One could also explore 
entanglement entropy in the presence of flavors for ABJM theory, see earlier results in 
[77, 78]. 

The entanglement entropy calculated using the ten-dimensional metric is in general quali¬ 
tatively different to that computed using the lower-dimensional Einstein metric. It is inter¬ 
esting to note however that the entanglement entropy computed using the top down metric 
for smeared solutions such as [77] seems to give answers which agree with E-theorem expec¬ 
tations. In such examples the ten-dimensional metric has a very special form, in which all 
warp factors depend only on the radial coordinate, and thus the lower-dimensional Einstein 
metric is simply related to the top down metric. It would be interesting to understand the 
relationship between the metrics in more detail, and to compare the entanglement entropy 
computed in [77] with what is obtained using the method developed here. 

The entanglement entropy is not the only interesting quantity which is computable from 
the effective lower-dimensional Einstein metric: correlation functions involving only the 
transverse traceless components of the held theory stress energy tensor can also be com¬ 
puted from perturbations of the Einstein metric. Kaluza-Klein holography allows such 
energy momentum tensor correlation functions to be accessed without computing the en¬ 
tire backreaction in ten dimensions. 

We were able to match the structure of all terms in the entanglement entropy for hnite mass 
quarks at zero density with held theory expectations, and the logarithmic divergences were 
matched exactly. Few held theory results exist for entanglement entropy at hnite tempera¬ 
ture and density; see the recent papers [79, 80] for discussions of specihc universal thermal 
corrections in conformal held theories. It would be interesting to explore whether conformal 
perturbation techniques analogous to those of [23] can be used to extract universal terms 
at hnite mass and density; just like the mass, the chemical potential breaks conformal 
invariance even at zero temperature but can be treated in conformal perturbation theory. 
Results on entanglement entropy using conformal perturbation theory in the context of 
higher spin theory can be found in [81, 82, 83]. 

In this work we have emphasised that the entanglement entropy and differential entropy 
are associated with the lower-dimensional Einstein metric, rather than the ten-dimensional 
metric, and that the latter can only be reconstructed given additional information. The 
emergence and reconstruction of the compact part of the bulk geometry is a longstanding 
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puzzle. If entanglement can indeed be used to reconstruct the bulk spacetime, then there 
must exist in the field theory a generalised measure of entanglement which captures the 
compact part of the geometry. Attempts to define such a quantity were made in [70, 73] 
and we will report elsewhere on extensions of these proposals. 
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A. Source terms in linear equations 


In this appendix we discuss the derivation of (5.19) from (5.17) using the results for the 
linearised field equations around an AdS^ x background given in [55]. The components 
of the Einstein equations in the non-compact directions are 
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Here and Vq denote covariant derivatives while and Dy denote the d’Alambertians; 
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Projecting this equation onto the zeroth spherical harmonic results in 
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where was defined in (5.7). 
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To obtain an equation for the Einstein metric perturbation we need to eliminate tt^ 
and The trace of the Einstein equation over the five sphere gives 

2 (□. - - m: + (A.3) 

Projecting this equation onto the zeroth spherical harmonic results in 

^(□,-32)7rO = Ac2orO. (A.4) 


The five-form self duality equation along the non-compact directions gives 
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1 


59[mCnpgr] ~ ^^^mnpqr daCbcde 2 ^^^ ^blg^emnpqr 
which projected onto the zeroth spherical harmonics gives 

rf) ffi — ^((hE\p _ 

'^‘^[m^npqr] ~ 'P 3 ^ )^mnpqr 

Inserting (A.4) and (A.6) into (A.2) then gives 

{Ce + 4)/l^„ = K?o(^mn + 


(A.5) 


(A.6) 


(A.7) 
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